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Transverse offset and angular misalignment (tilt) are serious causes of 
loss in multimode fiber splices. Our computation of these losses in multi- 
mode graded-index fibers reveals that the loss depends strongly on the 
power distribution in the fiber. We find that offsets of 0.1 of the core radius 
(or tilts of 0.1 of the fiber numerical aperture) cause a loss of 0.1 dB in 
the case of steady-state conditions, and between 0.34 and 0.38 dB if the 
power distribution is assumed as uniform. We compare these results with 
measured offset loss values and conclude that the steady-state distribution 
better reflects transmission line conditions than the uniform distribution. 

I. INTRODUCTION 

Fiber transmission lines will frequently have more than 10 splices 
between terminal points ; to estimate the overall line loss, we need 
a reliable figure for the splice loss to be expected under varying line 
conditions. The loss in multimode fiber splices depends on the power 
distribution among the modes in the fiber. The difficulty of correctly 
simulating line conditions in laboratory measurements or computa- 
tions has so far prevented the establishment of a reliable splice-loss 
estimate. Calculations have predicted as much as 0.4 dB for a 10- 
percent relative offset (offset divided by the core radius). These calcu- 
lations assumed an equal or uniform power distribution in all modes. 1 
Techniques used to measure the offset loss have attempted to simulate 
the equal-power condition by illuminating the full numerical aperture 
of the input fiber by a coherent source, 1 - 2 but the resulting loss values 
were between 0.1 and 0.2 dB for a 10-percent offset. The question has 
been raised whether some of the discrepancies are attributable to 
leaky rays. 

The following study considers lateral offset and angular misalign- 
ment (tilt) of splices made from perfectly similar fibers by joining 
them without air gap. As an introduction to the subject, we derive 
simplified offset-loss formulae for power-law profiles, considering both 
the absence and presence of leaky rays. We find that leaky rays are 
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most likely not the cause of the discrepancies. The main part of this 
study concerns the misalignment loss in square-law fibers when the 
mode-power distribution is nonuniform as a result of mode-conversion 
effects. Special attention is given to the steady-state power distribu- 
tion. 3 We find that the discrepancies mentioned earlier can be ex- 
plained in terms of differences in the relevant mode-power distributions. 

II. UNIFORM DISTRIBUTION 

To compute the offset loss, Miller 1 compares the number of rays 
that exit and enter at coinciding points of the two fiber end faces of 
an offset fiber joint (see Fig. la). The percentage of rays not accepted 
on entrance is the offset loss. As a first case, we assume all rays to be 
uniformly excited. Both fibers have the same dimensions and are of the 
same kind. In the case of step-index fibers, all rays traversing the 
overlapping core areas continue to propagate and the only loss occurs 
outside these areas. Offset graded-index fibers show a loss also in the 
overlapping areas where the local numerical aperture (lna) of the out- 
going fiber is smaller than the corresponding lna of the incoming 
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Fig. 1 — Areas of transmission in a fiber splice, (a) End-on view of the splice, (b) 
Area of overlap rearranged for computation of offset loss. 
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fiber.* This condition exists for example in the right-hand half of the 
area of overlap in Fig. la. In this area, the outgoing fiber accepts only 
those rays that fall within its lna; all others are lost. On the other 
hand — and this is important for our derivation — all rays in the out- 
going fiber are uniformly excited in this area since we have assumed 
that all rays of the incoming fiber were excited. 

In the left-hand half of the area of overlap, the incoming lna is 
smaller than the outgoing lna. Thus, all incoming rays are accepted 
there. On the other hand, the outgoing fiber is not fully excited. Since 
we assume that both fibers have the same characteristics, the number of 
incoming rays in the left half equals that of the accepted rays in the 
right half. If we arrange the two halves as in Fig. lb, we find that the 
number of rays lost equals the number of those rays that traverse the 
blank area of Fig. lb in a uniformly excited fiber. The width of this 
area is equal to the offset. 

Assume a fiber core of diameter 2a and define a radial coordinate r 
that is so normalized that r = 1 at the core periphery. Let n(r) be the 
core index and assume a constant index n(l) in the cladding. In the 
case of full and uniform excitation of all trapped modes, the rays 
traversing a point in the core are uniformly distributed over all angular 
directions within a cone that has the lna as its apex. The number of 
these rays is proportional to (lna) 2 = n 2 (r) — n 2 (l). The integration 
of this quantity over the blank area of Fig. lb, properly normalized, 
yields the offset loss. The integration can be performed analytically 
for step-index and parabolic profiles and leads to the results of Ref. 1. 

If the offset is a small fraction d of the fiber radius, we can approxi- 
mate the blank area by a rectangle of length 2a and width da. In addi- 
tion, we can ignore the variation of the lna across the width of that 
area. These assumptions simplify the integration so that the offset 
loss becomes simply 

, ( \j*(r) - n s (l)]dr 

L.--4 — (1) 

T / 2 (r) - n 2 {l)~\rdr 
Jo 

If the index follows a power law 

n*(r) = ng(l - 2Ar°) (2) 

in the core, the evaluation yields 

L, - ^«+f (3) 

jr a t" 1 



* The lna is proportional to the square root of the index difference at that point. 
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A comparison with the exact expressions for a = 1, 2, and °o shows 
that (3) represents the first term of an expansion of L t in powers of d. 
Remember that the result applies in the case of full and uniform excita- 
tion of all trapped rays, that is, those rays that are totally reflected 
inside the core. 

Leaky rays possess an evanescent cladding field in a region extending 
from the core periphery to a circle of radius A in the cladding. 4 The 
smaller A, the larger the radiative (leakage) loss 1(A). Stewart 4 
has found that the solid angle, which includes all rays having a loss 
smaller than 1(A), equals the solid angle comprising the trapped rays 
divided by {A 2 — r 2 )*. To obtain an upper bound to the influence of 
the leaky rays on the offset loss, we assume A = a in the following 
although, by doing so, we include leaky rays that have a significant 
leakage loss. 

Arguments analogous to the previous derivation lead to an offset 
loss 



L f = - 



f — 

A (i- 



r 2 )* 



dr 



L (1- 



(4) 



r 2 )* 



rdr 



for a uniform power distribution in both trapped and leaky modes. 
We obtain 



d 1-1/K 



with 



2 1 - K/(a + 1) 

x = 2- 2 ar 2 (f) / /r(«). 



(5) 
(6) 



Table I lists loss values for some exponents a. Figure 2 shows loss 
versus offset for a = 2. In that case, an offset of one tenth of the core 
radius causes a loss of 8.5 percent or 0.38 dB if only the trapped 
modes are considered and 7.5 percent or 0.34 dB if all leaky modes 
are taken into account as well. The discrepancy is small. 

Equations (1) and (4) are based on the assumption that at least 
one of the fibers of the joint is fully excited at any given point in the 
splice cross section. This assumption requires some scrutiny in the 
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Fig. 2— Offset loss vs relative offset for uniform power distribution according to 
eqs. (3) and (5). 

case of the leaky modes. In this case, the cone of rays traversing a 
given point in the core becomes elliptical, with the major axes oriented 
in radial and tangential direction. The offset juxtaposes cones with 
slightly different orientation. For (4) to hold, the larger of the cones 
must completely overlap the smaller cone. Without presenting details, 
we merely mention here that this condition is not satisfied, but that 
the resulting error is within the linear approximation in d. 

III. NONUNIFORM DISTRIBUTIONS 

Consider a ray that propagates in a fiber that has an index profile 
of the form (2). Assume that the ray forms an angle 9 with the fiber 
axis at a distance ra from the axis. The principal ray optics variables 
r and have no physical meaning in the wave optics picture; the 
combination 

(7) 



„ / , sin 2 0Y 



however, is directly related to the propagation constant of a fiber 
mode.* Given the power in each fiber mode, we can plot the power 
distribution p(R). Although the following computations can be per- 
formed for general power-law profiles, we restrict the derivation to 
the square-law profile in the interest of simplicity. The results for 
near square-law profiles are very similar. We introduce a normalized 



• ^J2 — (1 — p ! /k 2 nl)/2A, where /3 is the propagation constant and fc the vacuum 
wave number. R is the magnitude of a vector coordinate in phase space pointing to 
the locus of the mode. We call R, therefore, the mode coordinate. 
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angular coordinate 

p = sin 0/(2 A)*, (8) 

so that, for square-law profiles, 

jR 2 = r 2 + p 2 . (9) 

The fact that r and p are interchangeable in this relationship implies 
that offsets and tilts can be treated in the same way.* More specifi- 
cally, in case of a tilt angle 5 and a lateral offset d, we can introduce a 
single offset variable D so that 

D 2 = d 2 + sin 2 5/2A. (10) 

Before we calculate the offset loss, let us compute the total trans- 
mitted power P for a given power distribution p(R). We note from 
(9) that pdp = RdR and that p(R) =0 for r > R. Hence, 

P= C p{R)2icRdR f 2-icrdr (11) 

Jo Jo 



or 



P = 2ir* f p(R)R 3 dR. (12) 

Jo 

If the power distribution depends on the angle <f>, as indicated in Fig. 3, 
we obtain 

P = 4tt f R 3 dR f q(R,(j>) sin 2 <pd<j> (13) 

Jo Jo 

by introducing rdr = ydx with x = r cos 4> and y = r sin <p in (11). 

Now consider a lateral offset d. Assume the power distribution in 
the incoming fiber to be q(Ri) for R t < 1 and zero everywhere else. 
If we retain only linear terms of d as before, Ri = R — d cos <j>, where 
R is the mode coordinate of the outgoing fiber. The power distribution 
in terms of this coordinate becomes 

O raroosKS— l)/d] 

p(R) = - / q(R - d cos 4>) sin 2 <pd<f>. (14) 

7T Jo 

This expression is a modification of the second integral in (13), which 
takes into account that q = for > arcos £(R — l)/eQ. A similar 
derivation holds for a tilt angle 5 or a combination of offset and tilt. 
The result has the form (14) with d replaced by D of (10). It is con- 
venient to normalize q(R) for unit total power P; this is done in the 
following analyses. 



* The modes of square-law fibers occupy a phase-space volume of rotational 
symmetry. 
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Fig. 3 — Coordinate system used to compute offset loss. 



As a first example, assume a uniform input distribution 



2/V 2 




for 
for 



Ri < 1 
Ri > 1. 



Because of (14), the distribution in the splice becomes 

2 R- 1 |\ (R- I) 2 "!* 



Pi(R) = — 3 arcos — ^ — 



D 



b - VI 



(15) 



(16) 



Figure 4a shows q x and pi in the vicinity of R = 1 for D = 0.1. The 
offset loss is indicated by the cross-hatched area in Fig. 4a. It is 



f l+D 8D 

U = 2tt 2 / Pl (R)R 3 dR = 5- 

J l ov 



(17) 



in accordance with (3), which yields L t — 8d/37r for a = 2. 

Power transfer among the modes as a result of guide imperfections 
modifies the step distribution (15). If we assume uniform excitation 
at the input and model the power exchange by a diffusion process, 
we find that the step function (15) assumes the form 

1 -exp[(#,-- l)/0.17(z/z c )*] (18) 

a short distance z from the input. In this expression, z c is the coupling 
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Fig. 4 — Power distribution immediately ahead of the splice (q) and immediately 
behind the splice (p) for various power distributions plotted vs the mode coordinate 
R. The cross-hatched areas indicate splice loss; D = 0.1. 



length. 3 ' For our purposes, the approximation 



= 20+MJ (1 _^ )A 

[ 



92 = 



is satisfactory and valid as long as 

a = 0.17 (z/z e )l 



for 


Ri <1 - <T 


for 


1 - a < Ri < X 


for 


R<> 1 



(19) 



(20) 



is small compared to unity. Figures 4b and c show plots of 92 for 
a = 0.05 and a — 0.1 as well as the distribution pz(R) calculated from 
(14). An integration of the cross-hatched areas similar to (17) yields 
the loss L 2 . Although an analytical solution exists, we list here only the 
simpler approximate form 



9D 

U = — (1 + 2a) 

IT 



(f^Tfl+fs)" 1 . (") 



which is more useful for design consideration. Figure 5 shows Li and 
L 2 for a- = 0.05 and <r = 0.1 plotted vs the offset D. If, for example, 
z B = 200 m, a power distribution resembling that of Fig. 4b (a = 0.05) 



* The coupling length can be obtained by measuring the impulse response at a 
fiber length where the steady state has been reached. The coupling length is then equal 
to that length, multiplied by the square of the rms pulse width measured, and 
divided by the square of the rms width expected without mode conversion. 
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would exist 17.3 m from the input; similarly, a distribution resembling 
that of Fig. 4c would exist 34.6 m from the input (provided the input 
excitation is uniform). A 10-percent offset (or tilt) at these points 
would cause a loss of 0.14 dB or 0.09 dB respectively compared to 
0.38 dB if the distribution were uniform at the splice. The reason for 
the drastic decrease in offset loss can be found in the fact that the 
modes close to cutoff suffer the main part of the offset loss, but cease 
to carry a significant portion of the total power a short distance away 
from the (uniformly excited) input. 

For a > D, the last term in (21) dominates. The loss can then be 
written in the form 



L = |V(1W 



(22) 



where q'(l) is the derivative of g with respect to R at R = 1. In this 
case, the loss depends on the square of the offset and L(D) has zero 
slope at D — > 0. 

As z approaches z , the power distribution in the fiber resembles 
more and more the steady-state distribution. 3 Properly normalized, 
this distribution assumes the form 



q 3 = 



2tt 2 (w 2 - 4) 



J Q (uR)/Ji(u) 




for 
for 



R < 1, 
R > 1, 



(23) 



where J and J x are Bessel functions and u = 2.405 is the first root of 




RELATIVE OFFSET OR TILT 



Fig. 5 — Offset loss in dB vs relative offset or tilt for the four mode-power distribu- 
tions sketched in Fig. 4. 



OFFSET AND TILT LOSS 913 



Jo. Figure 4d shows g 3 as well as the distribution p 3 in the vicinity 
of R = 1. The fraction of p 3 that extends beyond R = 1 indicates the 
loss in the splice; it can be obtained from (22) by computing 93(1) 
with the help of (23). The resulting loss L 3 is also plotted in Fig. 5. 
The part of the power distribution p 3 that deviates from the steady- 
state g 3 in the region R < 1 is accepted by the outgoing fiber, but is 
dissipated later as a result of mode conversion until the steady state 
is reestablished. An integration of p 3 — q 3 analogous to (17) but ex- 
tending from R = 1 — D to R = I reveals that this loss contribution 
is equal to the loss suffered in the splice. Thus, the total loss can be 
obtained by inserting g 3 '(l) into (22), as mentioned earlier, and by 
doubling the result of this integration; the total loss caused by an 
offset under steady-state conditions is, therefore, 



L. = 



u A D 2 



8(u 2 - 4) 



(24) 



This loss is plotted in Fig. 6. A 10-percent offset or a tilt angle equal 
to one tenth of the fiber na causes 0.1 dB loss compared to 0.38 dB 
if the power is uniformly distributed in all trapped modes or 0.34 dB 
if leaky modes are present as well (Fig. 2). 




0.1 0.2 
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Fig. 6 — 'Steady-state loss vs relative offset. Points show loss measured for offset 
between two J-km fibers. 
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IV. COMPARISON WITH MEASUREMENTS 

Figure 6 also shows loss values measured in a splice connecting two 
2-km fibers. Details concerning the measurement are presented in 
Ref. 1. The coupling length characterizing these fibers was not mea- 
sured, but was probably of the order of the fiber length. The fact that 
the loss values are consistently lower than the predicted steady-state 
loss may indicate that the 5-km length was not quite sufficient for 
the steady-state power distribution to establish itself before and after 
the splice. 

Figure 7 shows loss values measured in a splice connecting two 3-m 
fiber pieces. 2 The fibers had an outside diameter of 125 /*m, a core 
diameter of 75 jum, and an na of 0.155. They were part of a cable struc- 
ture that may have introduced some mode conversion. The fiber input 
was excited by a HeNe laser. Focused on the core, the gaussian laser 
beam filled a numerical aperture of 0.25 (at 1/e 2 intensity). The solid 
line in Fig. 7 is a repetition of L 2 for o = 0.05 plotted in Fig. 5. Accord- 
ing to (20), this curve corresponds to the distribution present at a 
distance of 8.6 percent of the coupling length from the (uniformly 
excited) input. 

It is beyond the purpose of this paper to determine under what 
circumstances a coherent source produces a uniform power distribution 
at the fiber input. The results depicted in Fig. 7 suggest that the use 
of a laser beam for splice loss measurements is acceptable as long as the 
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Fig. 7 — 'Loss measured for offset between two 3-m fibers; solid line applies to 
conditions of Fig. 4b. 
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fiber na is adequately filled by the beam and a modest amount of mode 
conversion is present (or applied) between input and splice. Measure- 
ments under the conditions described yield loss values that are about 
50 percent too high at 10-percent offset and about equal to the steady 
state at 20-percent offset (or tilt). 

V. CONCLUSIONS 

We have attempted to show that the assumption of a uniform power 
distribution among all modes (or rays) is an idealization rarely ob- 
tained in fibers. The power distribution and particularly its behavior 
at the critical angle strongly affect the offset loss. If the offset loss 
calculation is based on the steady-state distribution, the loss at 10- 
percent relative offset (or tilt) is 3.8 times smaller than that predicted 
for a uniform distribution. Although the steady-state distribution does 
not necessarily exist at all splice points, we believe that it reflects 
transmission line conditions much better than a uniform distribution. 
We have compared these results with measurements by Miller 1 and 
Chinnock 2 and found that coherent excitation can provide offset loss 
values representative of line conditions if care is taken that the source 
fills the na of the fiber and a modest amount of mode conversion is 
effective ahead of the splice point. 
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